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Abstract. Results of identical pion correlations from the first year of data collection 
with the PHENIX detector at RHIC (7s^=130GeV) are presented. PHENIX has 
good particle identification using an electromagnetic calorimeter for timing, leading 
to identified pions from .2 to 1 GeV/c. This extends the range of previously measured 
correlation radii at this energy to (£;r)=633MeV/c. The beam energy dependence of 
the HBT radii are studied in depth and no significant dependence of the transverse radii 
is present. The longitudinal correlation length has a moderate energy dependence. 
Furthermore, theoretical predictions of R ou t/R.side severely underpredict the measured 
ratio, which is consistent with unity for all kj. The implications of these results are 
considered. 
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1. Theoretical Considerations 

1.1. The Physics of GGLP 

A distortion of the two particle probability distribution by the influence of Bose-Einstein 
statistics was first noted by the quartet of physicists Goldhaber, Goldhaber, Lee, and Pais 
in pp collisions in 1960 [ |]]. This observable was subsequently adopted in a variety of 
systems as a technique to probe the space-time history of the collision region. For a recent 
review see [ ||] . 

The normalized two particle probability distribution, P(pi,p2)/P(pi)P(pi), is given 
by the ratio of the two particle probability distribution to the square of the single particle 
distributions. This function is defined as the 'correlation function' and, for a static source 
with no final state interactions, is related to a Fourier transform of the source spatial distri- 
bution (p(r)) with respect to the momentum difference, q , C2 ~ 1 + p(q) 2 . A more general 
expression for the correlation function is: 
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C(ak)-l I \-fd^S( X , P )e-n 2 (1) 
C[q,k)-1 + lfd4xS ^ p)l2 (1) 

where k is the momentum average of the pair. 

Due to experimental statistical limits, it has to date been practically impossible to study 
the full six dimensional correlation function. Instead, in most cases we take advantage 
of the event averaged impact parameter to reduce our nominally azimuthally asymmetric 
source into a symmetric one []: (q,k) — > (q,kr,y) where kj is the component of k perpen- 
dicular to the beam. Historically, experimentalists have plotted C2 (q) while cutting on kj 
and y. 

1.2. A Note on Frames 

The appropriate choice of reference frame for the correlation function depends on our phys- 
ical assumptions^] of the source. Traditionally experimentalists and theorists in relativistic 
heavy ion collisions have chosen the longitudinally comoving system (LCMS) as the appro- 
priate frame and deconvolve q into components long, side and out. The 'long' component 
is that along the beam direction; the 'out' component is in the direction of the average 
transverse momentum of the pair; and the 'side' component is the complementary orthog- 
onal component, perpendicular to the transverse motion of the pair. This choice of frame 
assumes a longitudinally boost invariant source and has the very nice feature that a static 
long lived source results in R out ^> Rsuie- This feature is significant in searches for a quark 
gluon plasma where a phase transition is accompanied by a large latent heat and long lived 
source. 

However, if the source is expanding, the LCMS parametrization sees a Lorentz con- 
tracted outward radius (R^ t MS ), not the true coherence length of the source in its rest frame 
(i?™ )■ As such, measurements of the HBT radii are distorted not only by the fact that 
we measure only "lengths of homogeneity" but also by Lorentz boosts of the source in the 
outward direction. We can place some boundaries on the actual size of the outward radius 
in the source frame by noting that the largest value of R 0llt is measured in the pair center 
of mass frame (PCMS). All other frames which are related to the PCMS frame by a boost 
in the direction of the pair transverse momentum have a smaller measured outward radius, 
Rout MS /Jpr> where y pr is the gamma boost from the PCMS to the measured frame. A lower 
limit can also be placed on the true out radius if we are willing to assume that the source is 
not moving faster than the PCMS frame and is not going backward: R™^ E > R^ t MS . 

Consider a simple example of a source moving at (3 S = 0.7 toward the detector while 
emitting pairs with mean average momentum of (kj) = 300 MeV. In this case, if the out- 
ward radius measured in the LCMS frame is R L ™ S = 5 fm, then the radius measured in the 
PCMS is simply Y pr R L u t MS = 11.8 fm. Without a priori knowledge of the source velocity 
from, for example, assumptions from the singles spectra, we can only limit the true source 

1 see [ [3J] for a very interesting counter example. 

2 Much to our chagrin, mild assumption of the source are necessary even in the limit of infinite statistics due 
to the loss of information in Eq. [lj from the 7 dimensional source function (8 minus the mass constraint) to the 6 
dimensional correlation function. 
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PHENIX Detector - First Year Physics Run 




Fig. 1. The PHE^f^ experimefea'ftXh'fTguration in^rie first year of data taking, 
to be between 5 and 11.8 fm. However, the true source in its rest frame has a radius of 

2. Experimental Measures 

2.1. The PHENIX Experiment 

The PHENIX experiment has been described in detail elsewhere [ Q]. It is composed 
of 4 "arms": two arms perpendicular to the beam specializing in electron, photon and 
hadron identification ("central arms") and two arms for measuring muons ("muon arms"). 
A beams-eye view of the central arms is shown in Figure [j]. 

The central arms, covering a pseudo-rapidity region |r|| < 0.35, are composed of a 
variety of subdetectors used in concert for particle identification and momentum determi- 
nation. For the analysis in this note, we use detectors only on the west arms. The beam- 
beam counters (BBC) and zero degree calorimeters (ZDC) were used to trigger on inelastic 
interactions and determine the centrality of the collision; the drift chamber (DC), and pad 
chamber (PC) are used for three dimensional pattern recognition and momentum determi- 
nation; the association of the track with a cluster in the electromagnetic calorimeter (EMC) 
in conjunction with the BBC start time determines the velocity of the particle [ 

2.2. Analysis Details 

Centrality is defined by the measured zero degree neutral energy in the ZDC and the pro- 
duced particle multiplicity measured in the BBC [ ^j. For this analysis we consider events 
in the 30% most central range of the inelastic cross section. The mean centrality of all pairs 
in this sample is 10%. In the first year data sample this corresponds to 493K events after 
all offline cuts. 
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The mass of particles is determined in the analysis through a combination of momen- 
tum and velocity information. Momentum dependent bands are both measured experi- 
mentally and determined analytically from the momentum resolution of the DC and PC 
(8p/p = 0.6 ffi 3.6%/?) and timing resolution of the EMC-BBC pair (700 ps). Pions are 
defined as being within 1.5o of the pion mass peak and at least 2.5a from the kaon peak. 
The resulting dataset contains 3.1 million K + pairs and 3.3 million 7t~ pairs. The mixed 
background is constructed from the pair dataset through the random choice of two pions 
from pairs in different events. Due to the changing detector acceptance as a function of 
collision position we require that all mixed pairs come from events with a reconstructed 
BBC collision vertex within 1cm of each other. 

To further minimize the influence of detector induced artifacts in the correlation func- 
tion we introduce a number of cuts on both the real and mixed pair distributions. We require 
that all pairs are at least 2cm from each other in the drift chamber to suppress close-track 
inefficiencies and ghosts. Pairs that share the same EMC cluster are also removed from 
both the real and mixed samples. 

The correlation function is created by forming the ratio of real to mixed pairs versus the 
momentum difference of the pair. This raw correlation function is further corrected for the 
known Coulomb interaction between the two pions. The correction is performed by making 
the simplifying assumption of a Gaussian source in the pair center of mass frame and using 
an iterative procedure [ ^[ . The extracted radii have been studied as a function of both the 
radius (Rcoulomb) of the source and the strength of the coulomb interaction (Acaulomb) an d 
are found to vary by no more than .25 fm in any dimension versus reasonable variations 
of Rcoulomb (±20%) and Acoulomb (-5 — * 1-0). The results quoted below use the nominal 
extracted radius (Rcoulomb = Rim) an d full Coulomb correction (Acoulomb = 1-0)- 

Finally, the correlation function is fit using a log-likelihood MINUIT based minimiza- 
tion method to the functional form 

C 2 = l+A, Yl exp(-Ri 2 qi 2 ) 

i=x,y... 

where i= inv for the one dimensional fit, i = long, side, out for the LCMS fit, etc. 

Systematic errors enter primarily in the Coulomb correction and two track cuts. The 
total systematic errors of these effects is 8% for Ri on g and R s uie, and 4% for R out [ 

2.3. Results 

A projection of the %~ three dimensional correlation function is shown in Fig. [|with its 
corresponding fit. This particular plot includes the entire unbinned dataset. 

The kj dependence of the correlation function provides extra information on the dy- 
namics of the source so we divide our results into three bins in kj from 200-400 MeV, 
400-550 MeV and 550-1000 MeV. The corresponding (k T ) are 333 MeV, 472 MeV and 
633 MeV, respectively. In each kj bin we determine /?,„,,, as well as three dimensional fits 
in the LCMS and PCMS frame. The results are shown in Table [|. 
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Fig. 2. The three dimensional correlation function for Jt pairs versus qi ong , q S ide, and 

g'oujin both the LCMS frame (top) and pair center-of-mass frame (bottom). The data are 

plotted versus one momentum difference variable while requiring the other two to be less 

than 40 MeV/c. The lines correspond to the fit to the entire distribution. 

3. Interpretation 

The recent correlation results measured at RHIC, both by PHENIX [ |§ and STAR [|| 
have generated a great deal of interest and confusion in the theoretical community [ |ld| ]. 
Much of the non perturbative measurements at RHIC in the first year of data taking can be 
well described in the context of a hydrodynamic model. However, the correlation results 
have strained the parameters of these models to unphysical values, prompting a theoretical 
reevaluation. 



3.1. Comparisons 

Figure || places the measurements by PHENIX [ in the context of those at other energies 
and by other experiments [ ^, |llj [l2|, |l3| As is clear, the ^dependence of all three 
LCMS radii are very similar over a wide range of beam energies ( ^/snn =4.1 to 130 
GeV/c). Indeed the transverse radii are nearly identical over the entire energy range. The 
only clear energy dependence is in the longitudinal component which increases moderately 
as a function of energy. 

To quantify the differences in Ri on g as a function of beam energy, we fit each kj depen- 
dence to a hydrodynamically inspired equation [ |l5j [If]] f] : Ri ong = Aj ^Jmj. The results, 

overlaid with the data in Fig. f| are A = 3.32 ± 0.03, 2.9 ± 0.1 and 2.19 ± 0.05 fm-GeV^ 
for ^/s^=130, 17.3 and 4.9/4.1 GeV, respectively. 

3 We plot the ^Jsnn=A.9 [|l3|| and ^/.vjyjv=4. 1 [ |l4| l together due to their minimal energy difference. 
4 This functional form is supported by a theoretical approximation to first order of R/ ong versus T/mr- In this 
form, A = XqT where To is the proper hadronization time. 
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Table 1. The kj dependencies of the % + and % radii in the LCMS and PCMS frames. All 
momenta are in MeV and all radii are in fm. The errors are statistical only. 



k T (MeV) 

(k T ) 


200-400 
333 


400-550 
472 


550-1000 
633 


Rinv 

^■LCMS 

Rlong 

Rside 

Rout 
dPCMS 
"-out 


6.74 ±0.31 
0.423 ±0.037 
6.01 ±0.45 
4.81 ±0.30 
4.78 ±0.30 
11.35 ±0.69 


6.42 ± 0.46 
0.389 ±0.039 
4.76 ±0.35 
3.74 ±0.36 
3.76 ±0.26 
12.20 ±1.02 


3.46 ±0.46 
0.287 ±0.048 
2.97 ±0.38 
2.79 ±0.37 
2.59 ±0.46 
8.60±1.13 


Rinv 

K~ Rlong 

Rside 

Rout 
pPCMS 


6.00 ±0.30 
0.431 ±0.079 
5.69 ±0.76 
4.67 ±0.38 
4.69 ±0.58 
11.27 ±0.72 


5.96 ±0.41 
0.405 ±0.067 
4.77 ±0.49 
4.13 ±0.45 
3.75 ±0.40 
12.42± 1.18 


4.58 ±0.48 
0.353 ±0.062 
3.76±0.41 
3.22±0.35 
2.81 ±0.34 
11.89± 1.73 



There are at least two outstanding issues that have avoided easy physical descriptions: 
(1) the absolute magnitude of the radii and their shape in kj at RHIC is strikingly similar 
to lower energy measurements even though the particle multiplicity per unit rapidity in- 
creases by a factor of three over the same energy range; (2) the ratio of the transverse radii 
Rom /Rside is equal to 1 over all kj. 

The first puzzle, that of small radii, has been most troubling in the description of 
Rlong which has been irreproducible in hydrodynamic calculations, with predictions 20- 
50% higher than the data [ 0. 

The second puzzle, regarding the R ou t/ Rside ratio, is quantified in Figure ^ where we 
plot this ratio versus kj for both experiments at RHIC overlaid with theoretical predictions 
of a model that includes both a phase transition and a period of hadronic rescattering cal- 
culated with a microscopic transport code [ |7j. The model includes two scenarios with 
different critical temperatures and the authors note that a higher critical temperature leads 
to a longer period of rescattering in the hadronic phase. This in turn leads to a larger 
Rout /Rside ratio. The data are completely inconsistent with this theoretical prediction, and 
with any physically reasonable hydrodynamic theoretical calculation to date. [ [l(], |l8]| 

3.2. Hydrodynamical Fits 

The field of correlation measurements has adopted a method of interpretation based upon 
fits of the kj dependence of the radii. The effect of Bose-Einstein correlations on a boost- 
invariant, hydrodynamically expanding source can be calculated analytically [ [l9]] to first 
order 
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where R geom is the geometric radius of the source, m| = k\ + m|, P/ is the boost velocity 
(pV = p/ ^ p ) and T is the temperature. A fit of this functional form to the PHENIX data 

assuming T=125 MeV and pV = 0.69 [@, 0] Q is shown as a dashed line in Figure Q The 
extracted radius is R ge0 m =6.7±0.2 fm. 

Adding a second term in the expansion leads to the form [ |2j]] 

z? 2 

d 2 _ "geom 

^ -l +T1 2 ( l + ffi) (3) 

in which r|/ is the transverse rapidity boost T|r = T|/ r^ - ■ ^ e result °f a nt t0 me PHENIX 
data with this form yields R geom =8.1 ±0.3 fm. A fit of the STAR data to equation |] yields 
Rgeom =9.4±0.1 fm. These two fits are also shown in figure |i] All of these values for the 
geometric radius are much larger than the comparable ID rms radius of the Au nucleus. 

5 These values are taken from fits to the singles spectra for the 5%-15% centrality bin. 
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S^jjfrom identical pions for both 
R and PHENIX. The solid line is a 
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Taken together, the STAR and PHENIX data imply a steeper slope of R s Ue vs - 
then the fit, suggesting a higher p/T ratio. However, we should reserve our judgment 
of this issue until a full systematic study can be performed by both experiments with the 
larger second year data sets. At present, systematic errors in the experimental measures can 
completely account for the difference between the experiments. 

4. Conclusion 

We have presented the first year measurements of identical pion correlations as a function of 
kj by the PHENIX experiment at ^/saw=130 GeV. Comparisons to lower energies reveals a 
mild energy dependence in the longitudinal direction with no discernible transverse radius 
dependence. The ratio R ut/Rside is consistent with unity for all kj in contrast to theoretical 
predictions that include a phase transition. 



Acknowledgments 

The author would like to thank Dr. R.A. Soltz for enlightening discussions and unadulter- 
ated support. 

This work was performed under the auspices of the U.S. Department of Energy by 
the University of California, Lawrence Livermore National Laboratory under Contract No. 
W-7405-Eng-48. 



Particle Correlations with the PHENIX Experiment 



9 



References 

1. G. Goldhaber, S. Goldhaber, W. Lee, and A. Pais, Phys. Rev. 120, 200 (1960). 

2. U. Wiedemann and U. Heinz, Phys. Rep. 319, 145 (1999). 

3. M.A. Lisa, U. Heinz, and U.A. Wiedemann Phys. Lett B489 (2000) 287. 

4. D.P Morrison, Nucl. Phys. A638, 565c (1998), N. Saito, ibid, 575c (1998). 

5. J.T. Mitchell, et al, Nucl. Inst, and Meth. A482 (2002) 498. 

6. K. Adcox et al Phys. Rev. Lett. 86, 3500 (2001). 

7. M.D. Baker, Nucl. Phys. A 610, 213c (1996). 



8. K. Adcox et al, accepted for publication in Phys. Rev. Lett, nucl-ex/0201008 

9. C. Adler et al, Phys. Rev. Lett. 87, 082301 (2001). 

10. U. Heinz in these proceedings. 

11. I.G Beardenefa/., Phys. Rev. C58, 1656(1998). 

12. M.M. Aggarwal et al, Eur. Phys. Jour. C16, 445 (2000). 

13. R. Soltz, M.D. Baker, J.H. Lee, for the E802 Collaboration, Nucl. Phys. A661, 439 



(1999), L. Ahle et al, submitted to Phys. Rev. C (2002) |nucl-ex/0204001 . 

14. M. Lisa et al, Phys. Rev. Lett. 84, 2798 (2000). 

15. A.N. Makhlin and Y.M. Sinyukov, Z. Phys. C39, 69 (1988). 

16. U. Wiedemann, P. Scotto, U. Heinz, Phys. Rev. C 52, 918 (1996). 

17. S. Soff, S.A. Bass, and A. Dumitru, Phys. Rev. Lett. 86, 3981 (2001). 

18. D. Rischke and M. Gyulassy, Nucl. Phys. A608, 479 (1996). 

19. S. Chapman, J.R. Nix, U. Heinz, Phys. Rev. C52, 2694 (1995). 

20. J. Burward-Hoy, in these proceedings; K. Adcox et al, to be submitted to Phys. Rev. 
C(2002). 

21. U. Wiedemann, P. Scotto, U. Heinz, Phys. Rev. C53, 918 (1996). 



